Quantum Manifestation of Elastic Constants in Nanostructures 
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and outer surface with modified surface elastic constants 
llHl3j . This makes the overall mechanical properties of 
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Generally, there are two distinct effects in modifying the properties of low-dimensional nanostruc- 
tures: surface effect (SS) due to increased surface-volume ratio and quantum size effect (QSE) due 
to quantum confinement in reduced dimension. The SS has been widely shown to affect the elas- 
tic constants and mechanical properties of nanostructures. Here, using Pb nanofilm and graphene 
nanoribbon as model systems, we demonstrate the QSE on the elastic constants of nanostructures 
by first-principles calculations. We show that generally QSE is dominant in affecting the elastic 
constants of metallic nanostructures while SS is more pronounced in semiconductor and insulator 
nanostructures. Our findings have broad implications in quantum aspects of nanomechanics. 
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Mechanical properties of nanoscale solid structures 
are different from their bulk counterparts. It has been 
demonstrated experimentally that elastic moduli change 
their values as a function of the size of nanostruc- 
tures, such as the diameter of nanorods or thickness of 
nanoplates [l], Q ■ The general understanding is that such 
size dependence of elasticity has its physical origin in the 
elasticity variation at material's surface. It is well-known 
that surface has a different structure from underlying 
bulk due to bond breaking, surface relaxation and recon- 
struction [3H5|, which gives rise to excess surface energy 
and non-zero intrinsic surface stress • Consequently, 
the elastic constants of surface (which may include sev- 
eral atomic layers 0,01) are distinctively different from 
those of bulk. In a nanostructure, the surface-to- volume 
ratio continues to increase with the decreasing size, so 
that the overall elastic constants of the nanostructure 
will exhibit a strong size dependence. 

There have been many studies about the elastic con- 
stants of nanostructures focusing on the surface effect 
(SS). For example, experiments showed that Young's 
modulus of a thin film can either increase or decrease rela- 
tive to bulk when the film thickness approaches nanoscale 
d, 0] . Theoretically, it is found that the surface could 
decrease Young's modulus down to 2/3 of its bulk value 
from calculations usingharmonic or Lennard- Jones po- 
tential approximation [H [§[ . Another calculation found 
that Young's modulus of thin film varies as inverse of its 
thickness, which could go either larger or smaller than 
the bulk value, based on embedded atom method (EAM) 
and Stillinger- Weber potential [io| . EAM simulations 
of Cu film showed that Cu surface could become either 
stiffer or softer relative to bulk 



nanostructures distinctively different from those of their 
bulk counterpart. For example, the mechanical bending 
of nanofilms follows the modified Stoney [l2[ and Timo- 
shenko [13| formula rather than the classical formula for 
macroscopic thick films. 

Besides the SS, it is well known that there is another 
effect becomes increasingly prominent at the nanoscale 
to affect the properties of low-dimensional nanostruc- 
tures: quantum size effect (QSE) induced by quantum 
confinement. When the dimension of a nanostructure is 
reduced to be comparable to the electron Fermi wave- 
length, electrons become geometrically confined giving 
rise to quantized electronic states that change electronic 
energy, which in turn modify various properties of nanos- 
tructures by QSE, such as surface energy [14[, stabil- 



ity 15[ and magnetism f 16] . Therefore, it is reasonable 



11] . In general, elas- 



tic constants of nanostructures have been modeled by 
partitioning the structure into two parts of inner bulk 
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to expect the QSE to affect the mechanica l prope rties 
of nanostructures. A few recent theoretical [17H19| and 
experimental [20( studies have indeed shown the QSE 
causing quantum oscillations of surface (edge) stress in 
nanostructures. In general, however, despite the exten- 
sive study of the SS on the elastic constants of nanos- 
tructures little attention has been paid to 
the QSE. 

In this Letter, we demonstrate quantum manifestations 
of elastic constants in nanostructures induced by QSE 
using first-principles calculations. Using Pb nanofilms 
and graphene nanoribbons (GNRs) as model systems, 
we show that the Young's modulus and Poisson ratio 
of nanostructures can display an oscillatory dependence 
on size, i.e., the thickness of nanofilm and the width of 
nanoribbon. The main physical origin for such quan- 
tum oscillations of elastic constants is the QSE induced 
oscillation of electron density inside the nanostructure. 
Because electron Fermi wavelength is much shorter in a 
metal than in a semiconductor or insulator, generally the 
QSE dominates over the SS in affecting the elastic con- 
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FIG. 1: Schematic plot of computational supercell of (a) 
Pb(lll) film and (b) armchair GNR. Vacuum region was 
shown in outer blue box, and the thickness convention is 
shown in inner purple box. 




stants of metal nanostructures, while the reverse is true 
for semiconductor and insulator nanostructures. It is also 
important to point out that previous theoretical studies 
used empirical potentials which did not account 
for electronic effects. Consequently, the QSE on elastic 
constants will be missed in these earlier studies even if it 
were present, which calls for the need of first-principles 
methods. 

Our calculations are carried out using the density func- 
tional theory method as implemented in the VASP code 
[2l| with the projector augmented wave method [22[ and 
the Perdew-Burke- Ernzcrhof exchange-correlation func- 
tional 23]. As shown in Fig. 1, Pb(lll) film is modeled 
by a supercell slab set at the theoretical bulk lattice con- 
stant of 5.04A as the reference of strain-free state [24j . 
The slabs are separated by a vacuum thickness of 20A 
in z-direction, sampled by a 20x20x1 mesh in k-space. 
GNR is modeled by using similar super cell technique 
with a vacuum thickness of 20A in both y and z di- 
rections, sampled by a 10x1x1 mesh in k-space. All 
calculations used a plane- wave cutoff of 1.3 times of de- 
fault VASP value and the structure is optimized until the 
atomic forces converged to 1 meV/A. We extracted elas- 
tic constants of Yong's modulus (E) and Poisson ratio (u) 
from calculating stress-strain relations as a function of 
system size, for which we varied the Pb(lll) film thick- 
ness from 1 to 30 monolayers (MLs) and the armchair 
GNR (AGNR) width from 1 to 29 atomic rows. 

Young's modulus describes the stiffness of a material, 
defined as the ratio of tensile stress over tensile strain 



FIG. 2: (a) Young's modulus and (b) Poisson ratio of Pb(lll) 
film as a function of film thickness. The dashed lines show 
the fitted surface effect on Yong's modulus and Poisson ratio. 



(E = —). And Poisson ratio is defined as the ratio 
between biaxial transverse compressive strain over lon- 
gitudinal uniaxial tensile strain (y — Here, for 

Pb(lll) thin film, we apply biaxial compressive strain 
{s x ,y) in the film surface plane (normal to film surface), 
and calculate the strain induced film stress in the sur- 
face plane {a x ,y) (note that the intrinsic surface stress 
in the absence of strain is subtracted) and tensile strain 
(e z ) in the surface normal direction. Then we define the 



film Young's modulus as Et = 



Ao 



as uj 

for AGNR, we define E r = 



'-, as illustrated in Fig. 



and v r 



and Poisson ratio 
1(a). Similarly, 



'- , as illus- 



trated in Fig. 1(b). Another issue for systems with sur- 
face and edge is how to define their thickness and width. 
Here, we use the convention that the thickness (width) is 
set equal to the distance between the two outmost atomic 
planes (rows) plus one interlayer (inter-row) spacing, as 
show in Fig. 1(a) [Fig. 1(b)]. The interlayer spacing of 
Pb (111) film is 2.90A, and the inter-layer of armchair 
GNR is 3.35A. These same values are used throughout 
for consistency. 

Figure 2(a) and (b) shows the calculated Ef and Vf 
of Pb (111) film as a function film thickness from 1 to 
30 MLs, respectively. Clearly, both Ef and Vf show a 
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the elastic modulus of the film can be calculated as 



FIG. 3: Charge density distribution along z direction of 30Mb 
Pb(lll) film, (a) n m ax show the maximum charge den- 
sity within each atomic planes; (b) n m in show the minimum 
charge density in between atomic planes. 



strong odd-even oscillation with a beating pattern period 
of ~9MLs, manifesting the QSE. Overall, both oscillation 
patterns are very similar to those of surface energy (25j 
and surface stress of Pb (111) film. The oscillation 
amplitude for both Ef and vt decays slowly, remains 
to be strong even for ^30ML thick film. Apparently, 
they differ from commonly recognized trend that elas- 
tic constants change monotonically as a function of film 
thickness due to SS [il-fl3| . This clearly demonstrates the 
importance of QSE on modifying the elastic constants of 
nanofilms. 

To better understand the physical origin of quantum 
oscillations of elastic constants, we may consider a simple 
free electron gas model. The bulk modulus of a uniform 
electron gas of density n is [26( 



= ftW) 273 ^ K n 5/3 

3m 



(1) 



In a nanofilm, QSE modulates the electron density 
along the surface normal z-direction. Figure 3 shows the 
charge density distribution n{z) along the z-direction in 
the 30ML film. Clearly, both the maximum charge den- 
sity within an atomic plane (n m ax) and the minimum 
density in between two atomic planes (n m in) exhibit an 
odd-even periodic oscillation originated from the QSE 
along the z-direction of the film. Thus, approximately, 



If, , , 1 
Bf = — / B[n,z)dz oc — 



i(z) 5/3 dz 



(2) 



And the QSE modulated charge density distribution 
leads to the QSE modulated film modulus and similarly 
other elastic constants. 

Besides the QSE, the SS should be present also. If one 
considers the film has an elastic constants (Cb), thickness 
(d) and a surface layer of thickness (6) and surface elastic 
constants (C s ) the overall film elastic constant 

can be easily calculated as 



Cf - C b 



2{C S - C b )6 



(3) 



which shows an inverse linear dependence on the film 
thickness (~ 1/d) citell,12,13. Whether the surface be- 
comes harder or softer depends on C s . If C s > C&, C/ 
increases with decreasing d; if C s < Cb, Cf decreases 
with decreasing d. Now, if we pretend to ignore the QSE 
and use Eq. (3) to forcefully fit the calculated the Ef 
and Vf, we got the black dashed lines shown in Fig. 2, 
which in fact reflect the SS. From the fitting, we obtained 
Eh ^75 GPa and Vb ~0.5 at the bulk limit, which are in 
good agreement with experimental values of 80 GPa and 
i/=0.4 [27|. 

For comparison, we also preformed similar calculations 
for AGNRs. The reason for choosing AGNR is because 
it has been known that the AGNR exhibits interesting 
QSE effect in electron band structure 28j, with 1/3 be- 
ing metallic and 2/3 being semiconducting in a three- 
atomic-raw oscillation as a function of width, as well as 



similar oscillations in edge energy and edge stress [17 1 
(i.e., equivalent surface energy and stress in 2D). Figure 
4 shows the calculated Young's modulus (E r ) and Pois- 
son ratio (y r ) of AGNR as a function of width from 1 
to 29 atomic rows. Not too surprising, we see the sim- 
ilar quantum oscillations in both E r and v r induced by 
QSE, with a three-atomic-raw period as in edge energy 
and edge stress [17J . In addition, we also fit the data us- 
ing Eq. (3) to reveal the edge effect (i.e., the equivalent 
SS in 2D), shown as the dashed line in Fig. 4. From the 
fitting, we obtained the graphene Young's modulus and 
Poisson ratio to be 0.95 TPa and 0.16, in good agreement 
with previous works 2^, 30( | . 

It is interesting to compare the results of Pb(lll) 
film in Fig. 2 with those of AGNR in Fig. 4, to re- 
veal the relative importance of QSE versus SS. For Pb 
(111) film, the QSE modulation of elastic constant is so 
strong with an oscillation magnitude changing the elas- 
tic constants by ~100%, while the SS is less important 
changing the elastic constant by at most ~26%. In con- 
tract, for AGNR, the QSE induced oscillation magnitude 
is very small changing the Young's modulus by a maxi- 
mum of 2% and Poisson ratio by 11%, respectively, while 
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FIG. 4: (a) Young's modulus and (b) Poisson ratio of AGNR 
as a function of ribbon width. The dashed lines show the 
fitted surface effect on Young's modulus and Poisson ratio. 



the edge effect (the equivalent SS in 2D) changes them 
by as much as 7% and 340%, respectively. Therefore, 
we conclude that the QSE dominates over the SS in af- 
fecting the elastic constants of Pb nanofilms, while the 
SS dominates over the QSE in affecting the elastic con- 
stants of AGNRs. This can be generally understood as 
follows. First, consider a free electron gas model, the elec- 
tron Fermi wavelength (Xp) scales inversely with electron 
density (n) in a power- law (Xp oc n -1 / 3 for 3D electron 
gas and Xf oc iiT 1 !" 1 for 2D electron gas) [26], so that 
the Fermi wavelength is usually shorter in metals with a 
high electron density than in semiconductors and insula- 
tors with a low carrier density. Consequently, quantum 
confinement of electrons or carries, and hence the QSE is 
stronger in metal nanostructures than in semiconductor 
and insulator nanostructures. Second, the metal surfaces 
usually relax or reconstruct less than the semiconductor 
and insulator surfaces [3J, so that the SS is expected to 
be weaker in metal nanostructures than in semiconductor 
and insulator nanostructures. 

In summary, we have demonstrated that the QSE can 
have a profound effect in affecting the elastic constants 
of nanostructures, with interesting manifestation of size- 
dependent quantum oscillations, in addition to the mono- 
tonic size-dependence induced by SS that has been widely 



recognized before. Most importantly, we show that for 
metal nanostructures the QSE induced oscillations can 
be the most dominant effect to completely overwhelm the 
SS. Our findings shed important new light to our under- 
standing of the mechanical properties of nanostructures 
by adding interesting quantum aspects to nanomechanics 
with broad implications. 
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